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ABSTRACT 

The effects that interactions produce on gaiaxy disks and how they modify the subsequent 
formation of bars need to be distinguished to fully understand the relationship between bars 
and environment. To this aim we derive the bar fraction in three different environments ranging 
from the field to Virgo and Coma clusters, covering an unprecedentedly large range of galaxy 
luminosities (or, equivalently, stellar masses). We confirm that the fraction of barred galaxies 
strongly depends on galaxy luminosity. We also show that the difference between the bar fraction 
distributions as a function of galaxy luminosity (and mass) in the field and Coma cluster are 
statistically significant, with Virgo being an intermediate case. The fraction of barred galaxies 
shows a maximum of about 50% at M r ~ —20.5 in clusters, whereas the peak is shifted to 
M r ~ —19 in the field. We interpret this result as a variation of the effect of environment 
on bar formation depending on galaxy luminosity. We speculate that brighter disk galaxies are 
stable enough against interactions to keep their cold structure, thus, the interactions are able 
to trigger bar formation. For fainter galaxies the interactions become strong enough to heat up 
the disks inhibiting bar formation and even destroying the disks. Finally, we point out that the 
controversy regarding whether the bar fraction depends on environment could be resolved by 
taking into account the different luminosity ranges probed by the galaxy samples studied so far. 

Subject headings: galaxies: clusters: individual (Coma) — galaxies: clusters: individual (Virgo) - 
galaxies: evolution — galaxies: formation — galaxies: structure 



1. Introduction 

It is commonly accepted that galaxy bars spon- 
taneously form due to instabilities in dynamically 
cold disks (|Ostriker k. Peeblesl[l973l) . The growth 
rate of bars depends on the halo-to-disk mass ratio 
and the velocity dispersions of t he disk and halo 
( Athanassoula fc Sellwood 19861 ) and bars grow 
faster in massive and cold disks. Moreover, en- 
vironmental processes can regulate the life cycle 



of bars contributing both to their development by 
forcing disk inst abilities and to their destruction 
via disk heating (|Friedlilll999h . 

The role of environment in triggering the for- 
mation of bars has been a matter of discussion 
for a long time. A variety of methods to de- 
tect bars in galaxy disks and measure the local 
galaxy density have been adopted, but the results 
on the correlation between the bar fraction and 
environment are controversial. Since the early 
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work by iThompsonl (|198lh who claimed that the 
bar fraction of Coma galaxies increases toward 
the core of the cluster, similar res ults were found 
for the Virgo and Fo rnax clusters (Andersen 1996; 
Eskridge et al.l 2000) and for intermediate redshift 



clusters (jBarazza et al.ll2009l) . Furthermore, some 
authors meas ured an increase of the bar fraction in 
galax y pairs (iKumai et aTlll986UElmegreen et al.l 



Il990l: iGiuricin et al.l Il993t IVarela et al.l 12004 
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Recently, iBarwav et al. ( 2011 ) have found that 
fainter lenticular galaxies in clusters show a higher 
bar fraction than their counterparts in the field, 
and ISkibba et al.l (|2012l ). using clustering meth- 
ods, have argued that barred galaxies tend to 
populate high density environments. In contrast, 
other authors reported that environment does 
not play a major role in bar f ormation. Ac- 
cording to van den Berghl (2002 ), and later to 
Aguerri et al.l (|2009( ) and lLi et all (|2009l ). the bar 
fraction strongly depends on the properties of 
the host galax i es but not on their environment. 
Also Lee et all ( 2012 ) claimed that the bar frac- 
tion docs not depend on the environment when 
colour and central velocity dis persion are fixed. 
Mendez-Abreu et al.l (|2010 ) and lNair fc Abraham 
(l2010h jointed out that galaxy mass is the main 
parameter driving the bar formation and observed 
that bars a re hosted by galaxi e s in a tight range 
of masses. Martinez fe Muriell (|201lh found that 
the bar population does not significantly depend 
neither on group ma ss nor on the di s tance to the 
nearest neighbour. iGiordano et"aP (|201ll ) com- 
pared two carefully selected samples represen- 
tative of isolat e d and cluster galaxies, whereas 
Marinova et al.l ( 20121) investigated the bar frac- 
tion in lenticular galaxies across different envi- 
ronments which span two orders of magnitude in 
galaxy density. Neither of them found significant 
differences. 

Despite the fact that bars naturally appear in 
simulations of galaxy formation once a dynami- 
cally cold and rotational ly-supported disk is in 
place ( Athanassoulall200"5l ) . the use of bars as trac- 
ers of d isks has not been s o commonly adopted 
(but see iBarazza et al. I l2009l) . As for bars, the in- 
fluence of environment in the formation and evolu- 
tion of disks has been widely investigated. Phys- 
ical mechanisms taking place in galaxy clusters, 
such as harassment, ram pressure stripping, tidal 
effects and mergers, and starvation, were proposed 
for heavily transforming the morphology of disk s 
or even destroying them ( Boselli fc Gavazzl 2006'). 
Therefore, since the environment strongly affects 
the stellar and/or gaseous component of disks, it 
is customary to think that environmental effects 
have a direct influence on the formation of bars. 

In this Letter we derive the bar fraction in 
three different galaxy environments ranging from 
the field to the Virgo and Coma clusters. The 



unprecedentedly large range of luminosities (or, 
equivalently, stellar masses) covered by the dif- 
ferent galaxy samples we investigate allows us to 
distinguish the effects of environment in heating 
galaxy disks from those triggering bar formation. 
Throughout the paper we assume a fiat cosmology 
with J7m = 0.3, = 0.7, and a Hubble parameter 
H = 70 km s" 1 Mpc" 1 . 

2. Galaxy sample 

Four galaxy samples were selected in order to 
analyze three different galaxy environments: the 
field, the Virgo cluster, and the Coma cluster. 

The first galaxy sample (hereafter fieldl sam- 
ple) covers the low galaxy densi ty regime and was 
taken from Aguerri et al. (2009). They selected a 
volume-limited sample of ~ 2800 galaxies avail- 
able in the Slo an Digital Sky Survey Data Rele ase 
5 (SDSS-DR5. lAdelman-McCarthv et al.1l2007l) in 
the redshift range 0.01 < z < 0.04. The lo- 
cal projected number density around each sam- 
ple galaxy was computed using the distance of 
the galaxy to its fi fth nearest neighbour galaxy 
(|Balogh et al.ll2004l ). About 80% of the galax- 
ies are located in very low-density environments 
(£5 < 1 Mpc -2 ), whereas the remaining 20% 
reside in loose (£5 ~ 1 Mpc~ 2 ) and compact 
(£5 > 10 Mpc" 2 ) galaxy groups. Only the galax- 
ies in the lowest local density bin were included in 
the fieldl sample, thus containing 2389 galaxies in 
the magnitude range —24 < M r < —20. 

In order to extend our analysis to fainter field 
galax ies, we selected all the galaxies in the SDSS- 
DR7 (jAbazajian et al.ll2009l ) within 2500 < cz < 
3000 km s -1 and in the magnitude range —21 < 
M r < —13. This is our second galaxy sample 
(hereafter field2 sample) and it consists of 352 
galaxies. They are a subsample of the i s olated 
galaxies studied in ISanchez-Janssen et al. ( 201dh 
and represent a volume-limited sample complete 
out to M r = -15.5. 

The third galaxy sample (hereafter Virgo sam- 
ple) is composed of galaxies belonging to the Virgo 
cluster and was taken from Zarattini et al. (in 
prep.). Cluster members were selected to be ex- 
tended sources in the SDSS-DR7 with m r < 17.7, 
within one virial rad ius from the position of M87 
( Mamon et al.l 12004 ). and with recession veloci- 



ties available within ±2000 km s 1 with respect to 
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the redshift of Virgo cluster. This velocity range 
corresponds to a 3a cut on the vel ocity distribu- 



tion of the Virgo cluster galaxies (jBinggeli et al 



1987(). For galaxies without spectroscopy, the clus- 
ter members were selected to have a g — r color less 
tha n 0.2 mag above the value o f the red sequence 
fit (|Mendez-Abreu et all Eoiol ). Spurious back- 
ground objects were rejected by a visual inspection 
of their morphology f ollowi ng the prescriptions of 

chard fc Andreonl (|2008l ). The resulting Virgo 
sample is composed of 588 galaxies in the magni- 
tude range -22 < M r < -13. 

The fourth galaxy sample (hereafter Coma 
sample) includes galaxies be longing to the Coma 
cluste r and were selected by iMendez-Abreu et al. 
201dh. They are extended sourc es in SDSS-DR6 



(|Adelman-McCarthv etaH l2008h with m r < 21 
and within a 5-arcmin radius from the position 
of every pointing of the Hubble Space Telescope 
Advanced Camera for Su rveys (HST-ACS) C oma 
Cluster Treasury Survey (|Carter et al.ll2008l ). All 
the galaxies with spectroscopy have a recession 
velocity within ±3000 km s _1 with respect to the 
Coma cluster. This corresponds to a 3cr cut on the 
velocity distribution of the Coma cluster galaxies 
(jColless &: Dunnl ll996). As for the Virgo mem- 
bers, also the Coma photometric members were 
selected to have a g — r color less than 0.2 mag 
above the value of the red sequence fit and their 
morphology was inspected to reject background 
objects. The Coma sample consists of 169 galax- 
ies with —23 < M r < —14 and located mainly in 
the cluster center. 

In an effort to homogenize the dataset, the ax- 
ial ratios, g and r apparent magnitudes of all the 
galaxies in the four samples were retrieved from 
the latest available SDSS data release (SDSS III; 
lAihara et aD 1201 lh . For our analysis we consid- 
ered only the galaxies with b/a > 0.5 (a and b 
being the semimajor and semiminor axis lengths 
of the galaxies) in order to deal with projection 
effects. It is worth mentioning that the galaxy 
images of field2, Virgo, and Coma samples have 
a similar spatial resolution. Indeed, the farthest 
galaxies (i.e., those in the Coma sample) were an- 
alyzed by studying the HST images. At the dis- 
tance of the Coma cluster (100 Mpc), the resolu- 
tion of HST-ACS (O'.'l) corresponds to about 50 
pc. This gives essentially the same physical res- 
olution as SDSS observations have in the ficld2 



or in the Virgo cluster, and it allows us to re- 
solve bars down to sizes of fbar — 150 pc. The 
fieldl sample has a similar resolution for the clos- 
est galaxies but the mean resolution of th e sample 
is rbar — 1-3 kpc ( Aguerri et al. 20091) . How- 
ever, since the fieldl galaxies are the brightest 
(and, therefore, largest) galaxies in our sample, 
bars smaller than the resolution limit should be 
considered as nuclear bars and therefore they are 
not the subject of this paper. 

Another caveat could be due to the fact that 
on average cluster galaxies are redder than in the 
field. To avoid a possible color bias, we computed 
the stellar mass fr om the g — r color following 
Zibetti et all (|2009h . 



3. Identification of disks and detection of 
bars 

We adopted the morphological classification 
of the galaxies i n the SSDS-DR7 spectroscopi c 
sample given by Irluertas-Companv et all ( 201 lh ■ 
They divided the galaxies in four morphological 
classes (E, SO, Sab, Scd) based on t h e auto mated 



method of Irluertas-Companv et aL ( 2008 ) using 



learning machines to analyze the concentration 
and asymmetry parameters. The main new prop- 
erty of such a classification method is that a prob- 
ability of being in each of the four morphological 
classes is associated to each galaxy instead of as- 
signing it to a single class. We classified as a disk 
galaxy any galaxy with a probability of less than 
50% of being an elliptical. Due to the incom- 
pleteness of the SDSS spectroscopic sample and 
since several of our sample galaxies do not have 
any spectroscopic information, we remained with 
1604, 336, 228, and 44 disk galaxies in the fieldl, 
ficld2, Virgo, and Coma sample, respectively. 



Fo llowing our previous work s (|Aguerri et al 



20091 : IMendez-Abreu et al-llioToh . we detected the 
presence of bars in the different samples by visu- 
ally inspecting the galaxy images. We classified 
all the galaxies into strong barred, weakly barred, 
and unbarred. The classification of each sample 
was always performed by two of us. A caveat re- 
garding these criteria is that our distinction be- 
tween strong and weak bars is not directly related 
to the contribution of the bar to the total galaxy 
potential, but rather refers to a secure or possi- 
ble detection of a bar, respectively. Therefore, the 
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fraction of weak bars could be understood as an 
estimate of the uncertainty on the bar fraction. 

4. Bar fraction across different environ- 
ments 

For each galaxy sample we derived the ordinary 
bar fraction /d (calculated only for the disk galax- 
ies) and the overall bar fraction /t (calculated 
for all the galaxies independently of their Hubble 
type). Since bars can only be triggered in disks, /t 
combines the luminosity distribution of disk galax- 
ies with their probability of having a bar overcom- 
ing the problem of the identification of disk galax- 
ies. This is always a major concern in morphologi- 
cal classifications dealing with the measurement of 
bar fraction, /t allows us to probe a larger range 
of luminosities than /d . Figure Q] shows /d and 
/t as functions of the r-band absolute magnitude 
and mass of the galaxies in our four samples. To 
avoid issues related to the bin size, we applied a 
moving-average (boxcar) smoothing over the his- 
tograms using box widths of 1 mag and 0.5 dex and 
steps of 0.5 mag and 0.25 dex in magnitude and 
mass, respectively. The lower and upper bound- 
aries of the hatched areas correspond to the bar 
fraction calculated by considering only the strong 
(i.e., secure) bars and both the strong and weak 
(i.e., secure and uncertain) bars, respectively and 
including their statistical uncertainties. The latter 
were computed by estimating the confidence inter- 
vals on binomial po pulation proporti ons following 
the prescriptions by ICameron poll . The values 
of /d and /t for the samples of bright (fieldl) and 
faint field galaxies (field2) are in good agreement. 
Therefore, for studying the bar fraction distribu- 
tion we merged the two samples into a joint sam- 
ple of field galaxies (hereafter field sample) with 
magnitudes in the range —24 < M r < —15 (or 
10 7 < M*/M e < 10 12 ). 

We calculated the weighted mean, peak value, 
and corresponding errors in magnitude (and mass) 
of the bar fraction distributions of the samples by 
performing a series of 1000 Monte Carlo simula- 
tions taking into account the confidence intervals 
(Table [T]). We found that the bar fraction strongly 
depends on the galaxy lum inosity (or mass), as 
previo usly poi nted out by Mendez-Abreu et all 
(|2010l see also iNair fc Abraham! l2010h - In addi- 
tion, measuring /d and /t in different galaxy en- 



vironments allows us to show that the bar fraction 
distribution is different for galaxies living in the 
field or in clusters. According to the values given 
in TablcQ] this difference is statistically significant 
(> 68% confidence level) for the field and Coma 
samples, with Virgo sample being an intermediate 
case. 

The bar fraction of field galaxies peaks at M r ~ 
-19 (or M*/M Q ~ 10 9 ), whereas the largest 
fraction of barred galaxies in clusters is observed 
at M r ~ -20.5 (or M*/M & ~ 10 105 ). In 
both of them the maximum /d and /t are in 
the range 0.5 to 0.6. Moreover, there is marginal 
evidence in our cluster samples for two peaks in 
the overall bar fraction at M r ~ —20.5 and 
— 18.5, respective l y. Th is supports the claims of 
Nair fc Abraham! (|2010h for a bimodal distribu- 
tion of the bar fraction as a function of the lumi- 
nosity (or mass). 

For M r > -19.5 (or M*/M Q < 10 95 ), the 
bar fraction of the field galaxies is systematically 
larger than those of Virgo and Coma galaxies. We 
speculate about a hint that the bar fraction of 
the fainter galaxies in Virgo is also larger than 
that of Coma. On the contrary, the overall bar 
fraction is larger in cluster than in the field for 
-21.5 < M r < -19.5 (or 10 9 5 < M,/M < 
10 11 ). The bins of brighter (or more massive) 
galaxies of the cluster samples (M r < —21.5 or 
A4*/A40 > 10 11 ) suffer from incompleteness. 
This prevents a reliable comparison for both the 
ordinary and overall bar fractions with the field 
galaxies. 

5. Discussion and conclusions 

In order to study the influence of environment 
on bar formation, we derived the fraction of barred 
galaxies over a large range of galaxy luminosities 
(or, equivalently, stellar masses) in the two nearby 
benchmark clusters Virgo and Coma and com- 
pared it to that of a sample of field galaxies. Both 
the distributions of ordinary bar fraction /d and 
overall bar fraction /t were measured. The differ- 
ence of the bar fraction distributions as a function 
of galaxy luminosity (and mass) in the field and 
Coma cluster is found to be statistically significant 
(> 68% confidence level), with the Virgo cluster 
being an intermediate case. 

Since bars live in disks, these findings allow 
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Table 1: Weighted mean and peak value of the bar fraction distributions 



Bar Fraction Distribution 


Statistical Parameter 


Galaxy Property 


Field Sample 


Virgo Sample 


Coma Sample 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 



Overall 



Ordinary 



Mean 
Peak 
Moan 
Peak 



Luminosity 
Mass 

Luminosity 
Mass 

Luminosity 
Mass 

Luminosity 
Mass 



-19.45 ± 0.12 
9.61 ± 0.06(0.12) 

-19.07 ± 0.72 

9.49 ± 0.36(0.37) 
-18.81 ± 0.11 

9.50 ± 0.06(0.12) 

-19.07 ± 0.63 
9.50 ± 0.33(0.34) 



-19.50 ± 0.18 
9.94 ±0.07(0.12) 

-20.41 ± 0.81 
10.34 ± 0.44(0.45) 

-19.13 ± 0.25 
9.79 ± 0.13(0.16) 

-19.71 ± 0.67 
10.06 ± 0.38(0.39) 



-19.85 ± 0.25 
10.13 ± 0.09(0.13) 

-20.64 ± 0.45 
10.45 ± 0.24(0.26) 

-19.21 ± 0.29 
9.92 ± 0.12(0.16) 

-20.12 ± 0.78 
10.21 ± 0.37(0.38) 



Note. — uminosities and masses are given in r-band magnitudes and log (M* /Mo), res pectively. The erro rs on masses given 
in parentheses include the typical uncertainty in the optical mass-to-light ratios (0.1 dex. lZibetti et aUl2009l ). 



us to distinguish between the environmental pro- 
cesses inhibiting bar formation (heating) or even 
destroying the host disk from the processes trig- 
gering the disk instabilities which are responsible 
for bar formation. 

We interpret the decrease of /d and /t with 
decreasing galaxy luminosity observed for fainter 
(M r > —19 or, equivalently, less massive galax- 
ies {M^/Mq < 10 9 ) as due to the increase of 
the disk thickness. Indeed, Sanchez- Janssen et al 



(l2010h lave recently found that the minimum of 
the disk thickness distribution occurs at a char- 
acteristic mass M.*/ M.q ~ 2 x 10 9 (correspond- 
ing to M r ~ 18.5) below which low-mass galaxies 
start to be systematically thicker making it diffi- 
cult to develop a bar. We suggest that the val- 
ues of /d and /t in the field are systematically 
larger than those in Virgo and Coma because the 
low-mass galaxy disks in clusters are more eas- 
ily heated, or destroyed, by galaxy interactions 
and can not develop a bar. In the low luminosity 
regime nurture and nature are acting on galaxy 
disks in cluster and field, respectively. 

Since the values of /d and /t are larger for 
cluster galaxies with —21.5 < M r < —19.5 (or 
10 9 - 5 < M*/M Q < 10 11 ), we conclude that 
brighter disks are strong enough to survive inter- 
actions and form a bar. In the high luminosity 
regime nurture and nature are acting on bar for- 
mation in galaxy disks of cluster and field, respec- 
tively. 

From the theoretical point of view, environ- 
mental processes taking place in massive galaxy 



clusters do not seem to favor such an enhance- 
ment of the bar population. For instance, fast 
galaxy encounters are able to create bar-like fea- 
tures in galaxy disks depend ing mainly on the 
geom e try of the encounter iMastropietro et al 



2005; lAguerri fc Gonzalez-Garci'al l2009h . How- 



ever, these encounters are randomly oriented and 
can remove large amounts of stars from the galaxy 
outskirt s moving a galaxy from the bright to faint 
regime (jKormendv fc Benderl [2012h . Therefore, 
we suggest that the enhanced bar fraction ob- 
served in Virgo and Coma is not only due to pro- 
cesses taking place in clusters. 

One possible explanation could be related 
to the different morphological mixing in the 
field and clusters. In the luminosity range we 
probe, the galaxy populatio n in clusters is domi 
nated by l enticular galaxie s ( Binggeli et al. 19881) 



Masters et al 



Recen tly, ICameron et all (2010 ^ 
(|201ll ). and iMarinova et all (|2012[ ) have found 



that the fraction of barred galaxies is higher 
for lenticulars tha n for spirals. This is con trary 



to the result s of | Laurikainen et al.l (|2007l) and 
Aguerri et al. ( 20091) and also to the structural 
properties of lenticular galaxies (i.e., large bulge, 
thick disk, and little gas) which should inhibit the 
formation of bars. However, even if the higher 
fraction of barred galaxies in clusters was con- 
firmed to be due to lenticulars, this would point 
again to environmental processes acting in trigger- 
ing bar formation since the typical environmental 
mechanisms observed in clusters (i.e., harassment, 
high-speed encounters, or ram-pressure stripping) 
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Fig. 1. — Bar fraction distribution as function of the galaxy magnitudes (left panels) and masses (right 
panels). The bar fraction calculated using all the Hubble types (/t) and only the disk galaxies (/d) are 
plotted in the upper and bottom panels, respectively. The fieldl, field2, Virgo, and Coma samples are 
showed by the hatched orange, green, blue, and red areas, respectively. The grey arrow indicates the 
characteristic mass below w hich low-mass galaxies start to be systematically thicker (A4*/A4q ~ 2 x 10 9 ; 
Sanchez- Janssen et al. 20101) . 



should not be able to make disk galaxies, in gen- 
eral, more prone to bar instabilities. 

On the other hand, slow galaxy encounters have 
been suggested to be a very effic ient mechanism 
in inducing the formation of bars ( Noguchi 1988). 
These encounters are not likely to happen in rich 
clusters, like Virgo and Coma, but they are more 
common in galaxy groups. Therefore, we tenta- 
tively suggest that bar formation could preferen- 
tially occur in small-size galaxy groups before they 
fall into the cluster. If this is the case, we expect 
to observe a large bar fraction in galaxy groups. 
Moreover, some of the bars we are observing in 
clusters should represent the population of genuine 



bars induced by interactions and should have dif- 
ferent observable properties wit h respect to spon- 
taneou s bars as predicted by iMiwa fc Noguchil 
(jl998h . 



To summarize, we suggest a bimodality in the 
role of environment on bars, since interactions 
can trigger or inhibit the bar formation in galaxy 
disks depending on their luminosity. The disks of 
brighter galaxies are strong enough both to main- 
tain their cold structure and survive close inter- 
actions, which trigger the bar formation when the 
galaxies are probably in a pre-cluster stage. On 
the contrary, the disks of fainter galaxies are more 
fragile and hot. They can be heated or even de- 
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stroyed by galaxy interactions which inhibit the 
formation of bars. Classification methods (e.g., 
visual inspection, concentration cuts, color-based 
cuts) used to select a disk galaxy sample are known 
to add biases. These are usually considered as the 
main source of the different findings about the in- 
fluence of environment in bar formation. Our re- 
sults highlight the importance of studying galaxy 
samples which have been carefully selected in lu- 
minosity to avoid biases when dealing with bar 
statistics. We argue that most of the controver- 
sial results about the relationship between envi- 
ronment and bar fraction could be explained in 
terms of the different luminosity ranges covered 
by the galaxy samples studied so far. 
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